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XXV. On the Problem of Three Bodies. By the Rev. J. Challis, M.A., F.R.S., 
F.R.A.S., Plumian Professor of Astronomy and Experimental Philosophy in the 
University of Cambridge. 

Received May 15.— Read May 22, 1856. 

XHE determination of the motions of three bodies mutually attracting according 
to the law of gravity being a problem too complicated for exact solution, mathe- 
maticians have employed various methods of solving it approximately. It is well 
known that of these methods the one which appears to be the most obvious and 
direct, Introduces terms which may increase indefinitely with the time, and render 
the solution inapplicable to any observed case of motion. This difficulty occurs 
whether the problem be to find the perturbation of the moon's motion by the sun, or 
the perturbation of the motion of one planet by another, and the necessity of meeting 
or evading it has very much determined the courses which the solutions of these pro- 
blems have taken. In the theory of the moon's motion, Laplace, Pontecoulant, and 
others, have appealed to the results of observations of the motions of the moon's 
perigee and node, to justify the assumption of a form of solution which is not 
attended with the above-mentioned difficulty. Although this way of proceeding may 
lead to correct results, there can be no doubt that it is an abandonment of the prin- 
ciple of determining by analysis alone the form of development which is appropriate 
to the conditions of the problem. Again, in the theory of the motions of the planets, 
recourse is had on the same account to the method of the variation of parameters, 
more especially for determining the secular inequalities. Now it will perhaps be 
admitted that that method, elegant and exact though it be, is yet not indispensable, 
and that when it succeeds, there must be some direct method which would be equally 
successful and conduct to the same results. The discovery of such a method I have 
long considered to be a desideratum in the theory of gravitation, and having after 
much labour found out one by which the forms of the expressions for the radius-vector, 
longitude and latitude, and both the secular and the periodic inequalities, are evolved 
by the analysis alone, and which is applicable as well to the lunar as the planetary 
motions, I thought it might deserve the attention of the Royal Society. I propose 
in this communication to enter at length into the details of the method, and then to 
add a few remarks on its general principle, and to explain why, in common with the 
method of the variation of parameters, it succeeds in determining the motion of the 
apses of a disturbed orbit. 
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1. Let M represent the mass of the principal body, or its attraction at the unit of 
distance, m that of the body whose motion it is proposed to investigate, and m' that 
of the distm'bing body. The principal body is supposed to be at rest, and the rectan- 
gular coordinates and distances of the other two reckoned from its centre as origin, 
are respectively x, y, z, r and *', y', z', r', at the time t reckoned from a given epoch. 
Then if f* be put for M+m, and R for the expression 

m' {xaf + y^ + zJ) ^_ 

we have for determining the motion of m the known equations, 

"^+73+^=" ('•) 

d'^y ^y rfR 

^+^+-^=0 (2.) 

S+?+f =0 (3.) 

Analogous equations apply to the motion of m! as disturbed by m. The directions of 
the axes of coordinates are entirely arbitrary. Conceive, therefore, to be known at 
a given instant (To) the position of the plane passing through m in the direction of 
its motion at that instant, and through the centre of M, and let this plane be the 
plane oi xy. Conceive also to be known at a given instant (To) the position of the 
plane passing through m' in the direction of its motion at that instant, and through 
the centre of M, and let this plane make with the other the angle a. Also let the 
intersection of the two planes be the axis of x. 

2. Now since according to these suppositions, the body m would continue in the 
plane xy if the disturbing force of m' ceased at the time To, it is clear that the coordi- 
nate z at any time To4-'' is a small quantity of the order of the disturbing force. By 
multiplying the equations (1.), (2.), (3.) respectively by 2dx, 2dy, and 2dz, adding, 
and integrating, we obtain 

f +g+S-M^*+^=»> ■ • <^-> 

■Xp being the differential coefficient of R with respect to x, y, and z considered as 

functions of the time. But from what is said above, -^ and -j-'-jr are of the order 

of the square of the disturbing force. Hence as it is proposed to conduct the 
approximation according to the powers of the disturbing force, these terms in the 
first and second approximations must be omitted. Also, if 6 be the polar coordinate 

of m reckoned on the primitive plane of its orbit from the axis of x, and ^ &c. be 

neglected, x-=^r cos 6,y=r sin 6, r being now regarded as the projection of the distance 
on the plane oi xy. Consequently 
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{dR) _dR d^ dR dr . 
dt ~ d& 'dt'' dr 'dt ' 

and the equation (4.) becomes 

S+r.|-:-l^+2J(§.f+f .i;)..+C=0. ..... (5.) 

Again, the equations (1.) and (2.) give 

d^ d^ dR dR_ 

•^ dt^ ~y dt^ ~^^ dy~y dx—^- 

, ^d& 
d^y d'^x '^ dt , dR dR dR 

Hence, by integrating, ^^^~^~1^^^' ('^•) 

h being an arbitrary constant. Consequently, by substituting for -rr^ in (5.) from (6,), 

and neglecting the square of the disturbing force, 

dfi m 2f<. r^_'^Ji^dR ^C{dR di dR dr\^. 

dt^+r^-Y + ^—^J dd ^*~'^}\li'dt'^ dr 'dtj 

„^ CdR d& ,^ diCdR j^ Cd^9/CdR.\j, ,d&h 

»»* }WTf^*=It}T9'^*-]dr\}M^*r' and ^=^5 nearly. 

Hence it will be seen that to the first power of the disturbing force we have 

dr^ , h^ 2|x , ^ S{dH(QdR ,\ dR dr\ ,, ,^ , 

^+^-;+C=2j|^,(J^J^)-^.j,}rf^ (7.) 

The equations (6.) and (7.) are suitable for determining the forms of the develop- 
ments of r and & in terms of the time. 
The function R becomes, by neglecting utiz, 

iri{xo^ -\-yy^) nJ 

and reckoning & from the axis of x on the plane of the orbit of m' in its position at 
the time TJ,, we have to the same approximation, 

x'=r'cos^', y=r' sin ^' cos <y, s'=r'sin#sin<j/. 
Hence 

R=-j2 (cos^cos^'+sin^sin^cos») n^' 

'^ (r'2 + r^— 2?t' (cos fl cos 3' + sin fl sin 6' cos 0)))' 

If powers of a above the second be neglected, the following approximate value of R 

is obtained : 

■n mr ,^ ^. ml 

"^ ^ ' (r'2+r2-2?y'cos(6-6'))' 

— 2TO'sin i^sin ^sin^^l "Ta— "; ~r» r 

2\r'2 (r'^+r2-2n-'cos(fl-6')) / 
3 z2 
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3. The foregoing preliminaries having been gone through, the order in which the 
approximate integration is to be effected may now be stated. As the approximation 
is to proceed according to the disturbing force, the equations (6.) and (7.) must first 
be integrated omitting the terms involving R. We shall thus obtain values of r and 
6 as functions of ^ and arbitrary constants, just as in the case of the problem of two 
bodies, and these constants may be designated by the letters usually employed in that 
problem. As the exact values of the functions would be unsuitable for carrying on 
the approximation, they may be expanded in series proceeding according to the 
powers of the arbitrary constant e to as many terms as we please. In like manner 
the functions which express the values of r' and /in terms oft, maybe expanded 
according to the powers of e'. When these values of r, 6, r' and d' have been substi- 
tuted in the right-hand side of the equation (7.), that side becomes a function of t 
and constants ; and supposing the integrations indicated to have been effected, and 
the result to be Q, we shall have 

Q being a small quantity of the order of the disturbing force. 
Hence dt=dr-^-^,+^ -0+20^'^ 

Q^ &c. being neglected. In the second term we may substitute for r in terras of ^ 
from the first approximation, which gives 

,.(4'+l^_c)-=*(-J+^^-c)=£.*. 

Supposing, therefore, that by the first approximation ^.-a=/(#), we obtain 



dt(l+Qfit))= 



dr 



(-?+^c)* 



This equation being integrated, a relation is found between r, t and arbitrary con- 
stants, by means of which r is to be expressed in a series proceeding primarily 
according to the disturbing force, and subordinately according to the quantities e, e' 
and a>. This value of r is next to be substituted in the equation (6.), which, being 
put under the form 

shows that the right-band side then becomes a function of # and constants, &nd that 
by integration may be obtained in a series proceeding according to the same law of 
arrangement as the series for r. 
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The plane oi xy has hitherto been supposed to be coincident with the plane of »i's 
orbit at the time Tj. On this supposition values of r and 6 have been obtained, which 
fully take into account the first power of the disturbing force and the mutual incli- 
nation of the orbits. We are now at liberty to suppose the plane oi xy to have any 
other position making a small angle with the planes of the orbits of m and m' in their 
positions at the epochs T, and To, and the equation (3.), viz. 

d^s fuz dK 
^+78 + ^=0, 

may be employed for finding a series for % in terms of t. In the second term of this 
equation the value of r given by the second approximation is to be substituted, but 
in the third term it is only required to substitute the values of r, 6, r' and & given by 

the first approximation. Also for z and %' we may substitute in -j- the functions of 

t which express the values of these quantities on the supposition that the motions are 

ds 



undisturbed, and that they are referred to the new plane of xy. Thus -r- becomes 



a function of t and constants, and the above equation takes the form 

d^z 

which admits of being integrated only by successive approximations. 

The process by which it has been shown that r, 6 and z are approximated to, gives 
at the same time, mutatis mutandis, the values of r*, & and z' to the first power of the 
disturbing force of m. By means of these six quantities the approximation may be 
carried to terms inclusive of the square of the disturbing forces. 

Having thus exhibited the general scheme of this approximate solution of the 
Problem of Three Bodies, I proceed to exemplify its practicability. 

First Approximation. 

4. The first approximation, which omits the terms involving the disturbing force, 
and is therefore identical in form with the solution of the problem of two bodies, is 
obtained by integrating the equations 

dt=dr[-y^+-^-C'j , dd=^. 



V 



The first equation gives, by integration 



»^(<+T)=cos-^-i^/aV-(a-r)^ 

where for the sake of brevity a is put for t^, ^ for 1 §-, and n for _or-^. Let « 

be the constant introduced by the integration of the second equation, and in order 
to designate the constants in the present problem by the letters usually employed in 
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the elliptic theory, let s— ®- be put fornT. Then substituting p for w#+s— w, we have, 
as is known, the following expansions of r and 6 to the third power of e : 

-=1 — ecosp+2(l — cos2p)+-g-(cosj9— cos3/)) 

6=i-\-nt-\-2e sin/>+ "F sin 2p-\-g ( -5- sin 3j»— sin p ) . 

These values of r and 6 (excluding the terms involving e% for the sake of avoiding 
long calculations) will be employed in proceeding to the second approximation. It 
is evident that C^ h, ■ss and g must be regarded as the arbitrary constants of the inte- 
gration however far the approximation be carried, no other arbitrary quantities being 
introduced by the process. The quantities a and e are given functions of C and h, 
and at present they have no other signification. 

5. Before proceeding further, an inference may be drawn from the equation (7.) 
which will be useful hereafter. When the foregoing values of r and 6 are substituted 
in the right-hand side of that equation, the constant e will be a multiplier of that 
side, independently of any limitation of the orbit of w'. Now let, if possible, e=0. 

Then 1=— 2^, and the equation (7.) becomes 

Since the relation i/j^=h^C shows that C must be positive, it follows from the above 
equation that ^=0 and -=C, or that the orbit of m is a circle whose radius is equal 

to ^. But the orbit of m cannot be exactly a circle independently of the form and 

magnitude of the orbit of m', unless the disturbing force be indefinitely small. Con- 
sequently the supposition that e=0 draws with it the inference that the disturbing 
force vanishes. At the same time, the supposition that the disturbing force vanishes 
must leave e an arbitrary quantity, because on this supposition the problem is that of 
elliptic motion, and e is the eccentricity of the orbit. These conditions may be ana- 
lytically expressed by such an equation as e^=z^g-\-km', k being positive and of fixed 
value, and e^ being arbitrary*. 

Second Approximation. 

6. The first step towards expressing the right-hand side of the equation (7.) as a 
function of t, is to expand the quantity R in a series proceeding according to cosines 
of multiples of the arc ^—^. Let 

R=R„-}-R,cos(^— ^0+^2 cos 2(^-^')-f &c. 
=Ro+2.R.cos*(^— ^), 
the values of* being the integers 1, 2, 3, &c. Also let r=a{l-\-u), r'—a!{\-\-u'), 

* See Note (A) at the end of the paper. 
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d=:nt+i-j-v, 6'=n't-\-s'+v' ; and for the sake of brevity put q for nt+6—(n'i^-\-s'). 
Then m=— ecos/j+g— 2^08 2^, u'=—efcosp', 

6—&=q-\-v—v'=:q-\-2esmp—2e'sinp'. 

It is not necessary for our present purpose to employ expressions containing higher 

powers of e and e'. 

CdR 
7. It is next required to obtain an expression for | -jt <?^ in terms of t. Since each 

of the factors Ro, R„ R2, &c. is a function of r, r' and constants, it follows that 

-j^= — 2 • R,* sin s{6—d'). 
Also,if A5 represent the value of R,when a is substituted for r and a' for r', we have nearly 

' ' ' da da' 

. dK . ,dA, , 

' da '^ da! ^ 

Again, « sin «(^^^)=* sin «(§'+«— v') 

=s sin sq cos s{v—v')-\-s cos sq sin s{v—v'). 

Hence, by substituting the foregoing value oiv—v', and retaining only the first power 
of e and e', it will be found that 

s sin s{6—0')=s sin-y9'+e«^(sin {p+ sq)-\- sin (p—sq)^ 

—e's^{s,in(p'+sq)-{-sm{p'—sq)j. 

By supposing s to have all negative as well as positive integer values, this equation 
may be more briefly written thus : 

s sin s{d—ff)=-x sin «g+e*^sin (/'+*9') ~ eVsin {p'-{-sq). 

Now observing that ■$sin*5'C0S^=«sin (jo+*g'), because *=: + l, +2, +3, &c., the 
following result will be obtained by multiplying the foregoing values of R, and 

ss'ms{d—d'): 

R,ssms{d—^)=^-Y^insq+U'A,—^-^jesm{p+sq) 

Consequently, 

2^^ dt- -Js.R,* sin s{d-^)df 

. 2s^A,-as^ 

A, 



=^-^^ossq+X „+,(^_^,) ecos(y+.g) 

2^A, + a's-^ 
—2. — TTT rre' cos (p'+sq). 

n' + s{n—n') \r > i' 
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Hence, since ^=-—2e«* sin jB—5e^n* sin 2j9, the following equation will be obtained 
to terras of the second order with respect to e and e' : — 
^d%Cd^ ,, ^ 2n^As ... X 

+2-;iT^^) (2*'A -«*^^) e^ sin sq 

-2. («-T7|i:^(2*^A,-a.§) +^ A.) e^ sin i2p+sq) 

8. Similarly we have to express -t—^ as a function of t. Since -^ represents the 

partial differential coefficient of R with respect to r, 

dR dr _ dRs dr , , 

^■rf/=2--^-SCOS*(^-^), 

the values of* being 0, 1, 2, 3, &c. 

„ , „ . , rfA. ,dA, , ,.d^K «V .d^K , ,.d^K aV 



Hence 



dR,±_ du/dA, d^A, d^A, , A 

dr'dt dtKda'da^ "^ dadd )' 



Also TT=e« sin p+e*w sin 2j» nearly. 

Consequently to terms of the second order with respect to e. and e', 

dR dr dA, . , „ (dA, a d^A,\ . „ , , ^^^A, . , 

*^-^=«^W«'"^+«^''\^-2l^j sin2;,-ame« ^^sinpcosp'. 

Also to the first power of e and e', 

2cos*(^— ^)=2cos*g'— 2sin*5'.*(y— y') 

=2 cos *j+2e* (cos (jo+^g-) — cos (p—sq)) 
— 2eV (cos ip'-\-sq) — cos (^' — sq))) 

s being equal to 0, 1, 2. 3, &c. Or, if *=+0, +1, +2, &c. on the right-hand side of 
the equality, 

2 cos 5(^— ^) =cos sq-{-2es cos (p+sq) — 2e's cos (p'+sq). 

Hence, placing the terms corresponding to*=+0 apart from the others, and obser- 
ving that sin (p+sq) cosp' is equivalent to sin /? cos (p'+sq) when *= + 1, +2, &c., it 
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will be found that 

2 f .|=2am '^ ^np+ae^n (2 ^-a ^) sin 2^-2«eaV« '^, sinpcosp' 

+2.aen -^ sin {p-{-sq)— 2. a^ns-^ sin sq 
+2.(«e^«(l+*)^ ^ -^Jsm{2p+sq) 

— l.(2aee'ns -^+aa'ee'n -^A sin p cos (^ +*?), 

the values of* being now +1, +2, &c. 

9. We ai*e now prepared to express the right-hand side of the equation (7.) in 
terms of t. The results obtained in arts. 7 and 8 give, 

^ CdU d^ ^_ 

^ df}lA"-*~^ dr'dt~ 

-2am ^ sin;j-ae^« (2 %-a ^) sin 2p^2aea!^n J^, sin;? oos^ 

-2.esin(p+.?){;^A.+a«^j-2.e^sin(2p+.q).|^— ,A,+^^ 

-2.e^(sin(2p+*3)-sin*5).{^:j:^^(2*^A,- 

+2.ee' sin;, cos {p'+sq) .|^rp^z^(2.^A,+a'* ^') +2a«* ^+««'« ^ - 

For the sake of brevity of expression the following substitutions will be made : — 

_ 2n* . , dk, 

re— m' ' ' da 

^^ 5n,2 . , dK, a^n d^K 

M=;^:^A,+a«-^— y -^' 

After multiplying the right-hand side of the foregoing equation by dt, integrating, 
and substituting in the equation (7.), the following will be the result :— 

^ Lecos(j? + sg) , ^ Me^cos(2^+sg) . ^ /cos (2;? + ^g) _ cos sq \ 



Pee' 

;'+s(: 

MDCCCLVI. 4 A 
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Substituting Q for the sum of the terms on the right-hand side of this equation, and 
neglecting Q% &c., we have 

=*(-?+^^-C)-'-f(-5+^--C)-'Q 

because by the first approximation 

\ r^*^ r } d^ a*eVsin*5p(l+4ecosp) 

Consequently 

the constants a, e, s and m having the same signification as heretofore*. 

10. Before proceeding to eifect the integration above indicated, it will be right to 
remove certain analytical difficulties presented by the form of the equation. First, it 
may be urged that as Q contains the first power of e, the coefficient of the last terra 
might become infinite if e were indefinitely small, and the equation would no longer 
hold good. But it has already been proved (art. 5.) that e and the disturbing force 

vanish together, from which it follows that the quantity, - x disturbing force, may 

approach a finite value or zero as e diminishes. Again, it will be seen, if the quantity 
to be integrated be put under the form x(t)di, that the factor x(/) becomes infinite each 
time sin/?=0, and that the development fails for the values of t that satisfy this equa- 
tion. But it is well known that an analytical circumstance of this kind will not prevent 
our obtaining in the final analysis the correct development of r, provided the integration 
above indicated can be effected, and that the failure must admit of some interpretation 
relative to the proposed problem. Now it is not difficult to point out the significance 

of the failure in this instance. Let us suppose that cos"' ~-=p. Then the arc ?> 

can differ from /> only by a small quantity, and we have exactly r=a(l--ecos?i). 
Hence as the quantities a and e are absolutely constant, it would seem that the 
maximum and minimum values of r are a(l+e) and a(l— e) in every revolution of 
the disturbed body. This inference is manifestly untrue, and the reason that it has 
not been legitimately deduced is, that the above-mentioned failure occurs when r 
approaches a maximum or minimum value. The failure has, therefore, an important 
bearing on the problem, as showing that the maximum and minimum values of the 
radius- vector are not of constant magnitude. I proceed now with the integration. 

* See Note (B) at the end of the paper. 
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11. In order that our method of solution may be successful, the differential, 

-r-2-(l — 4ecos»), 

must admit of exact integration to terms inclusive of the second power of e. Now it 
will be seen by referring to the expression in art. 9 for which Q was substituted, that 
this integration depends on the following integrals, which are exact. The integra- 
tion can, therefore, be effected. 

/*cos p 1 4^ 

i;T;;2r(l— 4ecosM)c?<= z^-z-\ — cotp-\'Aet 

fcos2^ , coip Tn cos^ cosy + «' sin jB sin p' cosjb' 

fv cos(;) + gg) dt C^ / cos {p +sq) cos {p-sq) \ dt . „ . -^ 

J ^' n+s{n-n')' sm^p—^ ^'\n + s{n-n>yn-s{n-n')j%in''p l^—^> ^' ^^'J 

^-^-n^-^ln-r^r"^ [*=±1, ±2, &C.] 

r2- ^°I^f+t -^^<=f2.(g^^-^+ ^°^^f"t) ^^^ [.= 1,2,&C.] 

J n + s{n—n') sin^p J \n+s{n—n') ' n—s{n—n')Jsivrp ^ ' ' -> 

= -S- ».-,.(»_^y (^gff)+cosA-gcot;>) [«=±l,±2,&c.] 

r„ cos{2p+sq) _dt_C^ ( cos {2p + sq) cos {2p-sq) \ dt ^ o &o 1 

J ^'2w + s(n-n')'sin''i9— J '^\2n+s{n-n')'^2n-s{n-ri))siv?p L* — *> ^J ^'^'-J 

=-2-4^5::ii4r::^(?^+cos*?cotp) [*=+i, ±2, &c.] 

C^ / COS (2j3 + sq) _ cos sg \ dt _ ^ 2 cosQ^ + sg) 

}^'\2n + s{n—n') s{n-n')) aiiv'p~ s{n—n')[2n + s{n-n')y sin^ 

r / / , . ,v\ . . , , . \\ <^' V cos(»' + sff) 

j2.(«cosjpcos(y+*9)+(w'+*(«-«0)sini'sm(y+47))5j^=-2-— -j^^- 

Since the factors which multiply the two last differentials are different according 
as « is positive or negative, it is important to remark that these integrals have been 
obtained without giving to s the + and — signs. The factors which multiply the 
three preceding differentials do not contain s. 

From the results of these integrations the following equation is readily obtained : 

, /aedA,, „ , ^ , , aV d'A„ \ 1 

^ he cos sq 

n'^—s^{n—n')^' siap 

+2. -,3^^ (^(^^ Sin .y sin;,+cos .^ cosp)^ 
4 a2 
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2Me^ / 2n \ 1 

-'^- 4n^-s^{n-nlr [m^^^] '^" '^ S'H p+COS sq COSpj^ 

^ 2Ne^ cos jp+sq) 

' 's{n.—n'){2n+s{n-n')y sinp 

2ac^e^n d^A(, cosp' 
' n^—n'^' dada'' sinp 

.y Pe^ cos (j^ + sq) 

"^ n^—[n' + s{n.-n')y' sinp 

Since *=+!, ±2, &c., the following equalities are true: 

'!^ • cos sq cos p=^ 'S,.cos(p+sq) 

_ si n sq sinp -, cos (p +sq) 

2,'- - - r= — i. . 

s s 

Hence it will be found that the term in the above equation which involves cos (p+sq) 
is 

2 2 /' _ gLe^ I (M4-N)e° \ cos {p + sq) 

n—n' '\ s{n+s{n—n'))s{2n+s{n—n'))J sinp 

Consequently the equation (8.) becomes 

/ 3 dAo , 1 d^Ao\. , ,a-r 1 /-^-^ — ; - 

(«+;;^-^+2;i-&5-j<+s-«^=cos >— _-v«v-(«-r)^ 

-2^5^;7d^-t«-la(-2+7^cosp)+2;^.-^cos/)+^^::^.^^,cos 

^ L 1 2 « / 2L e (M + N)e \ 



Pe> 



+2-::5— 7:7-^^:7— :i7^. cos (/;'+j5)|. 



12. Before advancing to the next operation, our attention must be directed to the 
failure of the term containing the denominator w^— («'+*(«— «'))" in the case of 4=1. 
As the denominator vanishes for this value of s, it is necessary to retrace our steps 
and consider that case separately. Referring to the equation at the beginning of 
art. 9, it will be seen that if .y=l and P, represent the consequent value of P, the last 
term becomes P,ee'sin/»cos (/>'+?). Also, since 

we have \s\n p cos {p' -^ q)dt=\ smp cos {p-\-vf—xif')dt 

Thus the equiation (8.) will contain the term 

1 C^^ci^dtf 1 . , t \ 



2a^e^n, 
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which will give rise in the right-hand side of the equation at the end of art. 11, to 
the terms 

Hence, taking these terms into account, and putting the equation under the form 

^=cos{?.,+V»-C-?)+2-^?i^(U+eV+eW^^^ . . (9.) 
we shall have 

Ha dAo a^ d^Ao\ ?_ ^ / ~'^ , M+N \ 

^-\n da -^2n da^ )^^^P''^n-n' '^\sln + s{n-n')y s{2n + s{n-n'))r^^ '<P''^^9> 

It may be remarked, that in the terms containing the disturbing force we have put 
Pi for p and p[ for p', which is plainly allowable, because the reasoning might be 
repeated with these values in the place of p and p'. Also, it appears that the ex- 
pression for which W is substituted contains a term multiplied by t. This term 
might be included in p^, but it is more convenient to retain it in its present position. 
I proceed to develope r in terras of t by means of the equation (9.). 
13. This equation must give a result of this form, 



ae 



H and h representing respectively the terms which contain, and those which do not 
contain, the disturbing force. Hence, omitting A^ &c.. 

Consequently, putting g for the last term within the brackets of equation (9.), 
a—r f , /-; ^gs eUk ] 

— =cos|p,+ex/l-H^-;7Trp+4 

= cos (p,+es/T^W)+sm (p,-f-ev^T^;;^^)(^;;;=p-^) nearly. 
Hence li=cos{pi-\-e\/l—l¥) 

and h = sin {p,+e^/l — H.^) (--~=^—gj ■ 

By the first of these equations H may be developed in a series proceeding according 
to the powers of e, which will be found to be identical in form with the analogous 
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series in the elliptic theory. In the other equation terms involving e*x disturbing 
force are to be omitted, and cos j9, may therefore be put for H. Hence 

h=s\n (pi-\-e sin Pi){eh cot p^—g), 
or, to the same approximation as before, 

h=—g sinpi(l-\'2ecospi). 
Thus, since -=1— e(H+A)j we finally obtain 

^=l-eH+2ii(U+e(V+2Ucosp,)+e'W) (10.) 

It will be seen that in this process sinp^ has disappeared from the denominator. 

14. The expressions represented by U, V and W admit of simplifications which 
will render them more convenient for substitution in the equation (10.). In art. 9 
we have 



11 — «' ' ' da 



Hence ^.=--^--^Un-K^an{n-ri)^^ 

also, M=-^A,+a«^*-^^5^^=L+-i^A-^^?^'. 

Using these values, and putting hin—rl) under the form {s{n—ri)-\-n)—n, it will be 
readily found that 

2 /_ 2L M + N \ 

«-«' V s(n+s(w-»')) «(2« + s(n-w'))/ 

Again, since A, and -^ are homogeneous functions of a and a' of the dimensions — 1 
and —2 respectively, by a known theorem we have 

By these equations the differential coefiicients of A, with respect to a' may be elimi- 
nated. Thus by substituting in the expression for P in art. 9, and reducing, the 
following result will be obtained : — 

P=«rp^-{(2.^-*)«A,+((^-2.)(«-«')-«')«^j-«^'^^'- 
Hence, putting s=. 1, 

P.=2«A.-2„af'-a^«^'. 
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Also, since U= x^— S.-g — 37 Kscossq, 

' n da nr—s^\n—n'Y ^' 

_ -. Aae dKn ^ 2Le , , . 

2eUcos;),= -— -^ cos p,- 2 .^2_^a(^_^,)a cos {p,+sq). 

Consequently, by putting for H its known expression from the elliptic theory, sub- 

dA d^A d^A 

stituting —2 -^— « -^ for a'^^> and making use of the foregoing equalities, with 

the values of V and W given in art. 12, the equation (10.) becomes as follows : — 



1 dA, 



3e3, 






f 3e dAp e d'A^ ef 
~^\2an^ da'T'4n^ da^ "^BaV 



2n . , dA, 



(2kx—2a--^—a^ -^Vcos(Br— ?!r')+2w?sin(ra'— ©*)) Icosjp, 



1 ^ w— n' 



2^-*- n«-s>-n')^ ^°^*^ 



+^.2. 



2» . , dAs 

-,A,+a-r- 

^ — n da 



s\n-nf ^[n + s(n-n')y-n^ 



'{s^—s){n—n') — 3» 



dX.\ -"1 



n + s{n —k) 
3» 



/_2n . <?Aa 

\n-n<^''^"' da J 



+ 



w— w 



'> ^~ 2 ^a« 



.cos(p^+«9) 



dA. 



f(2s^-s)nA,+ ((s«-2s)(ii-n')-«')«^ , ,^^ 1 



'2 da"" 



In this equation * has the values +1, +2, +3, &c. in the terms containing cosAg 
and cos (p,H-«?), and the values +0,-1, +2, &c. in the term containing cosip^-^-sq). 

On comparing this expression for the radius-vector with that obtained by Laplace*, 
terms will be found in the latter identical with all the above, with the exception of 
those that contain e'cos;?,; and there are other terms to which none of the above cor- 
respond. These are only differences in form, arising from difference in the processes 
of integration. It is chiefly important to remark, that in the foregoing expression 
for r there is no term containing ent as a factor. The signification of that which 
contains ^nt will be presently considered. 

16. Having obtained the development of the radius-vector, it is easy to infer that 
of the longitude {6) from the equation (6.), viz. 

dt~r^~r^J di "*' 

and from the value o{\-m<^^ ^^ terms of t, which has been already found. Putting 
r,+hr for r, and taking h- to represent the terms multiplied by the disturbing force, 

* M6canique C6L part 1. liv. ii. No. 50, 
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we shall have to the same approximation as before, 

=;3-|3(l+3ecosjo,)Sr-^(l+2ecosp,lJ-^(/#; 

The development of the term S-^dt will be of the same form as in the elliptic theory. 
An equation obtained in art. 7 gives, 

1 CdR ,^ 1 ^ As , 1 -c' ' da r i \ 



^ 28^K^a's^ 
+2^«2- rJ+s{n-/) e'c08{p,^-sq); 

J 2slA, + «'.^* 

Since the relations of the constants h, a, n and e are expressed in our probleni in the 
same manner as in the elliptic theory, we have h=-nc^\/\ — e^. Hence we may put 
«a^ for A in the terms involving the disturbing force. Consequently, omitting e% &c., 

2m a 1 dh-s 
J^§recos^,^^=J^(-^-^-^2- ^._^.(^_^,). cossqyco^p^dt 

Also, to the same approximation, 

^{^^ ^^)^ cosjp,^#=Jj 2-^ cos sg cos />.,«f^ 

=;;i2-7 ,.^ / rNrSin(»,+,yo). 

In all these equations the values of* are +1, +;2, +3, &c. 

Putting now, for the sake of brevity, D for 2A, — 2a^-a''^, and F, G and tt for 
the coefficients of cos sq, cos {p,-^sq), cos {p',+sq) respectively in the expression for ~, 

a 

and observing that 
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the following equation is readily obtained : 

+2^sin (;?;+?) 

wF . , 2neG • / , v , 2ne'H • / / , x 

~ ':^JZ 1\ Sin sq-\—r? ; r^ Sin {p.+sq)-\ — tt-i ', Tx SHI (» +*9), 

a's{n-vl) 2 a2(»i + s(w— n')) ^' ^ o^(n.' + s(n— »')) 

the value *=1 being excluded from the last terra. 

Substituting the results of these integrations, the following will be found to be the 
value of ^: 

6=^+ [n+--^y+^ Sin 2p,-\-i {j sin 3;.,-sin p,") 

+P+^^ "^-2^« -^+4^4cos(z.-.!r')-2«f sm(^-..'))|sin/?, 

_i_J_v / 2"^ A. \ . 

~'"2«2 ^' Vs(n-«')~s(«-«')V^'" *^ 



— sS'-TT-T a1 — 7-, — ; r-^H-2nHlsin {'p',+sq). 



where * has the values +0, —1, +2, +3, &c. in the last term, and the values +1, 

+2, +3, &c. in the other terms, and A_j has the same value as A,. 

16. The expressions for r and 6 obtained in arts. 14 and 15 may be put under 

forms more compact, and more convenient for drawing inferences, by making the 

following substitutions : 

1 dAo ae^ 

A—a— ^--^+-2- 

n=^-(^- ^ + 2;i-^+4^COs(«r-,.)j^ 

N=«+- ~, SO that p,=Nt+e-n 

' na da ' ^' 

^'^=-2i(^>-«^+^)- 
MDCCCLVI. 4 B 
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Thus we shall have, 

r=K-a{'E.-ef) cos (N^+s-D)-^ cos 2(N#+8-n)+&c. 

-^2.Fcos*(N<+s-N7-0 

+^ 2 . G cos {*(N^f e-N'#-g') +N/+ s- n} 

+^2.IIcos{(a— l)(N#+g-N'#-s')+N#+s-n'}. 

And putting F', G', H' for the coefficients of sin*^, sin {p,-\-sq) and sin {'p,+sq) in 
the development of 6, we have 

^=g+N#4-2(E+e/")sin(N#+g-n)+^sin2(N#+e-n)+&c. 
+e'^sin(N#+g-n')4-2^2.F'sin*(N#+e-NV-0 



+-2 2.G'sin{*(N;+s-N'#-0+N#+s-n} 



+^2.H'sin{(*-l)(N;+s-NV-g')+N;+s-n'}. 

17. I proceed next to draw some conclusions from these values of the radius- 
vector and longitude. 

(1) The quantity A is the non-periodic part of the radius-vector, and being equal 

to a — 2 ■^+Y ^^ ^ function of given quantities and arbitrary constants. A is, 

therefore, invariable. It may also be remarked, that as the value of r may be put to 
the same approximation under the form 

/^('~A'^j('~T+^'^'- + P^'"'^^'^* terms j, 

1 dK 
the quantity «— ^ -j^ is approximately the mean distance. Thus, so far as this 

approximation shows, the mean distance is invariable. 

(2) The mean motion is necessarily the factor of the non-periodic term N^ in the 
development of 6. Hence 

Mean motion =N=«-J — -r-^- 

' na da 

For the reason just adduced, the mean motion is thus proved to be invariable. 

As the two quantities A and N are functions of a and e, they are by consequence 
functions of the arbitrary constants h and C. Hence, if the values of the non- 
periodic part of the radius-vector and the mean motion be deduced from observa- 
tion, the constants a and e, or h and C, become known. 

(3) The quantity z, being simply an arbitrary constant, is invariable. Analogous 
considerations apply to the mean distance, mean motion, and the epoch («') of the 
orbit of m! as disturbed by m. 



PROFESSOR CHALLIS ON THE PROBLEM OF THREE BODIES. 541 

(4) The expressions for E and H show that these quantities contain terms which 
have t for a factor, and may therefore increase indefinitely. This circumstance creates 
no difficulty with regard to 11, which is always part of a circular arc affected by a 
sine or cosine. But as E appears as a coefficient, it might seem that the develop- 
ments of r and & contain terms which admit of indefinite increase. It must, how- 
ever, be observed, that according to the remark made at the end of art. 12, the 
function that has given rise to these terms is really affected by a cosine, and that 
they have their origin in the development of that function in terms arranged accord- 
ing to the disturbing force. 

The following considerations will enable us to obtain, at least approximately, the 
periodic functions of which 11 and E are partial developments. Whatever functions 
the complete values of 11 and E are of t, they may be expanded in series of the form 
cs-)-j3#-f-yf-|-&c., the two first terms of which are already determined. Hence 

^^=/3-t-27/+&c. 5 = /3'-h2y7+&c. 

Let t be indefinitely small. Then substituting 

na da ' 2n dc? ' 
we shall have strictly the values of -^ and -^ for the epoch at which t commences, 

VIZ- ^=l^-4^cos(z^~t.) ■^=-^^,sm(t.-t^')- 

Now if t commenced at a different epoch, we should obtain for -^ and -rr the same 

expressions as those above, but different in value, because by hypothesis these 
differential coefficients vary with the time. The changes of value, which in actual 
cases take place very slowly, are due to changes in the eccentricities, and in the lon- 
gitudes of the apses, of the two orbits, and will be very approximately taken into 
account by substituting in the above equations for e, e', w and w', the variable quan- 

titles E, E', n and 11'. Like considerations apply to the values of -^ and -^. Thus 

we shall have four differential equations, by the simultaneous integration of which the 
four quantities may be obtained as periodic functions of the time. The arbitrary 
constants introduced by the integration are determined by the known values of the 
functions when #=0. These periodic functions are to be substituted for E, E', 11 
and 11', wherever these quantities occur in the developments of r, r', & and ^; The 
four equations just mentioned are identical with those obtained by the method of the 
variation of parameters for determining the eccentricities and longitudes of the 
apses. It is worthy of remark, that in both methods the changes of the eccentricities 
and of the longitudes of the apses which are due to the disturbances, are taken into 
account in calculating the changes themselves, so that the approximation does in 

4 B 2 
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fact extend beyond the first power of the disturbing force, so far as it relates to these 
two elements. 

If the approximation be made to include generally the square of the disturbing 
force, and the values of r and 6 in art. 16, and the like values of r* and 6', be used 
for that purpose, terms may arise containing coefficients which have f for a factor. 
These terms may be converted into periodic functions of the time by the application 
of the principles exhibited above, but in that case the differential equations by which 
E, E', n and U' are found will be of the second order, and the periodic functions 
will be more completely determined. 

The inferences (I), (2), (3) and (4) respecting the secular variations of the elements, 
although obtained in a manner quite new, agree exactly with those deduced from 
previous solutions of the same problem. 

18. Having now obtained the developments of r, 6, r* and 6', inclusive of both 
periodic and secular inequalities, to an extent which is suflScient for most of the appli- 
cations of the Planetary Theory, I shall reserve for a future opportunity the investi- 
gation of the inequalities in latitude, and shall then take occasion to show in detail 
how this method adapts itself to the determination of the motions of the moon. At 
present I propose, in concluding this memoir, to make a few general remarks on the 
Problem of Three Bodies. 

It has been already observed, that the solution here adopted introduces no terms 
containing ent in the coefl5cients. These terms are to be distinguished from those 
whose coefficients contain e'nt, which, as we have seen, have reference to seculai- 
variations of the eccentricity and of the motion of the apse, and would vanish with 
the eccentricity of the orbit of the disturbing body. The former relate to the motion 
itself of the apse, and are not peculiar to the Problem of Three Bodies, occurring in 
fact in cases where the force is directed to a fixed centre. To illustrate this remark, 

let us suppose the force directed to a fixed centre to be -^—(ju'r. Then, the differen- 
tial equation for finding the orbit being 

let this equation be integrated by successive approximations, first neglecting the last 
tei-m, and then substituting in that term the value of r given by the first approxima- 
tion. By this process a term containing t in the coefficient will be introduced, and 
the motion of the apse will fail of being ascertained. But if, instead of this process, 
the equation dr'^ ^^ 2w 

be obtained, ?ind its approximate integration be conducted according to the powers 
of |!*', no such term will arise, and the motion of the apse will be determined. The 
latter process is-exactly analogous to steps employed in the foregoing solution of the 
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Problem of Three Bodies. The difference of the analytical results of the two methods 

may be thus explained. The equation obtained by putting ;^ = 0, viz. A^— 2/-(<r+Cr° 

— l*'r*=0, may be shown to have three positive roots if C be positive, so that analyti- 
cally there are three apsidal distances. The first method, by embracing the third 
apsidai distance (no step being taken to exclude it), applies to the other two only in 
an expanded form, the expansion giving rise to terms containing the factor t. The 
other, by commencing the approximation with the equation 

^+72 7 + c=o, 

restricts the application of the solution to the part of the curve which has two apsidal 
distances, and accordingly finds the function of t which in the other method is ex- 
panded. The method of the variation of parameters, by the very nature of the pro- 
cess, restricts the analysis in the Problem of Three Bodies to two apsidal distances, 
and on this account is successful in determining the motion of the apse. 

Again, I think it important to remark that the solution of the Problem of Three 
Bodies, as here proposed, applies equally to the Lunar and the Planetary Theories. 
The Problem of the Moon's motion does not differ in the analytical treatment it 
requires, from that of the motion of a Planet. In the one case as well as the other the 
approximation ought to be conducted primarily according to the disturbing force, 
which is assumed to be small compared to the principal force, and secondarily 
according to the form of the orbit, which is assumed to differ little from a circle. It 
is not necessary to take account of the ratio of «' to n in arranging the developments, 
but only in estimating the magnitude and importance of the terms resulting from 
the integrations. The possibility of effecting the integrations is the proper proof of 
the correctness of the process, and of its being adequate to give the development which 
is alone appropriate to the question, and which must result from every process that 
is in all respects legitimate. After making any assumption respecting the analytical 
form of the solution (as in the Lunar Theory is done by introducing the constants c 
and g), there can be no certainty that the solution will not at some stage become 
empirical. Probably the reason that the process which succeeds for a planet has not 
been applied to the moon, is the difficulty of extending it to the square and higher 
powers of the disturbance (which in the Lunar Theory it is necessary to take into 
account), and of embracing in the same operation both the periodic and the secular 
inequalities. The method I have exhibited in this communication appears to meet 
this difficulty by evolving simultaneously both kinds of inequalities by a process 
which obviously may be extended to higher powers of the eccentricity and the dis- 
turbing force. Such an extension would require nothing more than great labour in 
executing the analytical details. 
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NOTE (A). 

Received July 7, 1856. 

It has been shown in art. 5 from a priori considerations, that if the constant e=0, 
the disturbing force vanishes, and if the disturbing force vanishes, e is arbitrary. 
Hence it appears from the development of the radius-vector in art. 16, that the eccen- 
tricity of the disturbed orbit and the disturbing force are related in such a manner 
that if the eccentricity =0, the disturbing force vanishes, and if the disturbing 
force =0, the eccentricity remains arbitrary. The particular relation which satisfies 
these conditions ought plainly to result from the solution of the Problem of Three 
Bodies, and it may, therefore, be worth while to inquire how far such a result can 
be deduced from the integrations effected in the foregoing approximate solution. 
Now the expressions for the variations of the eccentricity and of the longitude of 
the apse obtained in art. 17, are identical with those given by the method of the 
variation of parameters. Hence for the present purpose I may make use of the 
deductions from those expressions which are usually given in treatises on the Plane- 
tary Theory. Referring to Pratt's ' Mechanical Philosophy,' art. 385, we have the 
equations, 

Dg =— (2BD-J-CD') . . . (1.) EA=— (2BE4-CE') . . . . (2.) 

D'g=— (2B'D'-f CD) . . . (3.) E'A=^(2B'E'4-CE) . . . (4.) 

e;=D=' -}- £="-!- 2DE cos {(^-/i)#-{- A-/} (5.) 

g^_Df2_l_Ef2_|_2D'E'cos{(g-A);-fft-/} (6.) 

n«»n'B+«'a'mB' 

In these equations B, B' and C are known quantities independent of the eccentricities 
and longitudes of the apses, e, and e' are respectively the eccentricities of the orbits 
of tiie disturbed and disturbing bodies, k and / are arbitrary quantities, and D, D', E, E' 
are also arbitrary, excepting so far as they are connected by the first four equations. 
Let, if possible, e,=0 independently of the time. Then it follows from (5.) that D=0 
and E=0. Hence, since e\ does not consequently vanish, it appears by (6.) that D' 
and E' do not on this supposition both vanish, and, therefoi-e, by the first or second 
equation, that ?w'=0. Again, let m'=0. Then by (1.) and (2.), D^=0 and Eh=0, 
and by (7.), one of the quantities g and h vanishes. Hence one of the quantities D 
and E vanishes and the other remains arbitrary. Hence also e, is arbitrary. These 
results confirm the reasoning in art. 6. 



g ov h= +- {{nam'B—n'a'ma'y-{-nnaa'mm'C^}'' (7-) 
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NOTE (B). 

The following method of obtaining an expression for dt equivalent to that in 
art. 9, was communicated to me by Sir John Lubbock after the reading of my Paper, 
and appears to be well worthy of being recorded in connection with the process of 
solution I have adopted, as, on a resumption of the reasoning for the purpose of 
carrying the approximation farther, it might considerably abbreviate the analytical 
details. 

'• If dt=-\^-rdv, and v be taken for the independent variable, the equation 



becomes 
and if 






ar 



r= a — ae cos v + cos vXQ. sin vdv — sin v I Q cos vdv, 

iff- ^~7' COS V I Sin V L I 

and v^cos^M- 1- iQsinvd'y — \Q cos vdv}- 

u , /a J /a rdr + r sin vdv fO, sin vdv + r cos vdvfQ, cos vdv 

V ft *^ f* {a^e^—a—r+{cosvJ'Q,sinvdv~sinvJ'Q,cosvdvY}'^ 

Neglecting powers of Q above the first, 

, /ft r rdr rdv {sin vfO, sin vdv + cos vJ'Q cos vdv) 

~"^ /* l(aV-(«-r)2)* («V-(a-r)2)* 

r{a— r)dr {cos vfQ, sin vdv — sin vfOi cos vdv)l 
{a^e^-{a-rY)i I 

=\ /- ■ I 1-\ : — (sin v\ Q sin vdv+cos v \ Q cos vdv) 

7* COS vdv m * 

-|-— ^"2— (coswjQsinfcfo— sin v\Q.co^vdv) 

='\ / - n tH ^~2~ I Q sin vdv U 

which equation is true to all powers of the eccentricities and inclinations, v being the 
eccentric anomaly." 



